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Abstract 

The average persistent current (I) of diffusive electrons in the Hartree-Fock 
approximation is derived in a simple non-diagrammatic picture. The Fourier 
expansion directly reflects the winding number decomposition of the diffusive 
motion around the ring. One recovers the results of Ambegaokar and Eckern, 
and Schmid. Moreover one finds an expression for (I) which is valid beyond 
the diffusive regime. 
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The physics of persistent currents in mesoscopic isolated rings pierced by an Aharonov- 
Bohm flux <j) , has attracted a lot of interest on the description of the thermodynamic proper- 
ties of mesoscopic metals, both in the non-interacting picture or in the presence of electron- 
electron interactions. The description of the interactions is a complicated task: although 
several attempts have been made to describe the role of the interactions in one-dimensional 
or few-channel rings, with the help of analytical arguments or numerical calculations [|TJ, 
the diffusive nature of the electronic motion, which is probably essential in the experiments 
on metallic rings with finite thickness ||-f§], has been treated only in two series of papers 
by Ambegaokar and Eckern and Schmid |5|||. They calculated the average persistent cur- 
rent in the Hartree-Fock apprximation where the interaction is treated perturbatively in a 
diagrammatic picture. Although this calculation provides an average current smaller that 
experimentally observed ||, one can believe that it contains the essential physical ingredi- 
ent namely weakly interacting diffusive electrons. Thus it may be interesting to simplify as 
much as possible the calculation in order to possibly generalize it, or to compare it with nu- 
merical calculations. In the above papers, the average persistent current is calculated with a 
diagrammatic technique where the diffusive motion is described by a Cooperon pole with a 
wave vector quantized by periodic boundary conditions. Then, by Poisson summation, this 
sum over diffusive modes is transformed into a sum over harmonics with periodicity 0o/2 
where (f> is the flux quantum. 

In this short note, we propose a simple derivation where the average current is directly 
related to the return probability for a diffusive particle. Then this return probability is 
expanded according to the winding number of the diffusive motion around the ring, which 
gives directly access to the harmonics expansion of the current. Moreover, we obtain a 
general expression which can be used beyond the diffusive regime. 

Consider a quasi-one dimensional ring of perimeter L and of transverse section S, in 
which the motion is supposed to be diffusive along the perimeter and uniform along the 
transverse section. 

The first step is to write the total energy Et in the Hartree-Fock approximation: 
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E T = E° T + J2j t/(r-r')|^(r')ri^(r)rrfrdr' (1) 
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where is the total energy in the absence of interaction. In the lowest order in the interac- 
tion parameter, the states ipi are the states of the non-interacting system. The summation 
is over filled energy levels, cxj is the spin of a state ipi- 
Considering that the Coulomb interaction U(r — r') is screened in the metallic regime, 
it is replaced by U(r — r') = U5(r — r') where U = ^T\e 2 jq\ F and qTF is the Thomas- Fermi 
wave vector. Replacing the interaction in eq.[l| by a 5 function is certainly correct as long 
as the Thomas- Fermi wave length is smaller than the mean free path l e : qrple ^> 1- The 
interaction being now considered as 5-like, it is quite easy to see that the Fock term has the 
same structure as the Hartree term. Introducing the local density n(r) = J2i iV'iMI 2 ) the 
total energy can be rewritten: 

E T = E° T + U J n 2 {r)dv - - J n 2 (r)dv (2) 

The Fock term is half the Hartree contribution because of the constraint on the spin and its 
sign is opposite because of exchange of particles. The local density n(r) can be expressed in 
terms of the Green function: n(r) = (—1/tt) Jq f ImG R (r, r,u)dou so that the average current 
is given by 

(Ie-M)) = ir~) (3) 
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where G R (G A ) is the retarded (advanced) Green function. The product (G R (r,r)G A (r,r)) 
simply expresses the probability to go from some point r to itself ||. More precisely, for a 
particle at energy E 

P(E, u) = ^(G R (v, r, E + uj/2)G A (r, r, E - u/2)) (4) 

is the Fourier transform of the return probability P(E,t) after a time t: 



P(E,t) = — [ P{E,u)e~ iuJt duj (5) 
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Because of disorder average, this return probability is independent of the position r. The 
current can now be expressed directly in terms of P(E, t). Neglecting the energy dependence 
(P(E,t) = P(t)), one gets: 

= jj-fM* (6 ) 

Q = LS is the volume. Since the relation @ is exact, the expression (||) is quite general 
and is valid beyond the diffusive regime. In the classical approximation for the diffusive 
regime (l e <C L), the diffusion probability is the solution of a classical diffusion equation 
DAP = dP/dt where D is the diffusion coefficient taken here at the Fermi level, D = f^Z e /3. 
It is given by 

p,( T r ' t) = - e -|r-r'|2/4Dt (J) 

In the geometry of a quasi-one-dimensional ring, the return probability can thus be expanded 
according to the winding number of the diffusive motion: 

P (*' ® = o /7-n7 E t 1 + cos(47rm^)] (8) 

where p = 4>/4>o- The second term, of importance here, results from the phase interference 
between time-reversed paths in the semi-classical approximation, each path accumulating a 
phase ±2nmip In zero flux, the return probability is twice the classical one. This 

expansion according to the winding number gives directly the Fourier decomposition of the 
current: 

(J e _ e (0)) = ^2l m sin(47rmy?) (9) 

m 

with 
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We have introduced the Thouless energy E c = hD/L 2 . Defining a dimensionless winding 

number w by w 2 = m 2 /4E c t, I m can be simply rewritten as 

64 Up E c [°° 2 o 
I m = —/= — 2 w e dw 



= 16^^ (11) 

00 

which is the result of Ambegaokar and Eckern || and Schmid ||. 

We finish with the calculation of the current at finite temperature T where two Fermi 
factors /(e) have to be introduced in eq.|. Doing the standard substitution / f(e)g(e)de = 
2m Y^unOfan) where u n = (2n + l)vrT, the average current at finite T is straightforwardly 
given by: 

(I e ^ e ) = -QUp -^-4n 2 T 2 P{iu n -i<) 

TX Off) , 

=- nu *lhlj£m p(tA)it (12) 

By introducing the same dimensionless winding number w as above, the harmonics J m 
are given by 

w^w4f*' v ^>"sk (13) 
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where T m = E c /m 2 is the effective temperature associated with the winding number m 
and 9 m = T/T m . Although the integrand is quite different, this temperature dependence is 
identical to the one found by Ambegaokar and Eckern ||. It directly expresses the current 
in terms of a temperature square average of a winding number. 

In conclusion, we have calculated the average persistent current in the first order of the 
Hartree-Fock approximation. By writing the current directly in terms of the winding number 
decomposition of the return probability, we have avoided the use of the diagrammatic calcu- 
lation and directly found the harmonic expansion of the current. Although this calculation 
is reminiscent of the semiclassical description of the spectral correlations and of the average 
current of non- interacting particles in the canonical ensemble []T2| , it does not use any semi- 



classical sum rule, since here the correlation function of interest can be exactly written in 



terms of the return probability without any approximation. Moreover, we have written an 
expression for the average current (eqs|||12|) which is valid beyond the diffusive regime and 
may be used for example in the ballistic regime where interesting magnetic response have 
also been observed . 
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